Solving Nonlinear Parabolic Equations by a Strongly Implicit Finite-Difference Scheme

We discuss the numerical solution of nonlinear parabolic partial differential equations, exhibiting finite speed of propagation, via a strongly implicit finite-difference scheme with formal truncation error $\mathcal{O}\left[(\Delta x)^2 + (\Delta t)^2 \right]$. Our application of interest is the spreading of viscous gravity currents in the study of which these type of differential equations arise. Viscous gravity currents are low Reynolds number (viscous forces dominate inertial forces) flow phenomena in which a dense, viscous fluid displaces a lighter (usually immiscible) fluid. The fluids may be confined by the sidewalls of a channel or propagate in an unconfined two-dimensional (or axisymmetric three-dimensional) geometry. Under the lubrication approximation, the mathematical description of the spreading of these fluids reduces to solving the so-called thin-film equation for the current's shape $h(x,t)$. To solve such nonlinear parabolic equations we propose a finite-difference scheme based on the Crank--Nicolson idea. We implement the scheme for problems involving a single spatial coordinate (i.e., two-dimensional, axisymmetric or spherically-symmetric three-dimensional currents) on an equispaced but staggered grid. We benchmark the scheme against analytical solutions and highlight its strong numerical stability by specifically considering the spreading of non-Newtonian power-law fluids in a variable-width confined channel-like geometry (a "Hele-Shaw cell") subject to a given mass conservation/balance constraint. We show that this constraint can be implemented by re-expressing it as nonlinear flux boundary conditions on the domain's endpoints. Then, we show numerically that the scheme achieves its full second-order accuracy in space and time. We also highlight through numerical simulations how the proposed scheme accurately respects the mass conservation/balance constraint.Comment: 36 pages, 9 figures, Springer book class; v2 includes improvements and corrections; to appear as a contribution in "Applied Wave Mathematics II

Nonlinearity and Topology

The interplay of nonlinearity and topology results in many novel and emergent properties across a number of physical systems such as chiral magnets, nematic liquid crystals, Bose-Einstein condensates, photonics, high energy physics, etc. It also results in a wide variety of topological defects such as solitons, vortices, skyrmions, merons, hopfions, monopoles to name just a few. Interaction among and collision of these nontrivial defects itself is a topic of great interest. Curvature and underlying geometry also affect the shape, interaction and behavior of these defects. Such properties can be studied using techniques such as, e.g. the Bogomolnyi decomposition. Some applications of this interplay, e.g. in nonreciprocal photonics as well as topological materials such as Dirac and Weyl semimetals, are also elucidated

The IceCube Neutrino Observatory - Contributions to ICRC 2017 Part VI: IceCube-Gen2, the Next Generation Neutrino Observatory

Papers on research & development towards IceCube-Gen2, the next generation neutrino observatory at South Pole, submitted to the 35th International Cosmic Ray Conference (ICRC 2017, Busan, South Korea) by the IceCube-Gen2 Collaboration

Multi-messenger observations of a binary neutron star merger

On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ~1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40+8-8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 Mo. An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ~40 Mpc) less than 11 hours after the merger by the One- Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ~10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ~9 and ~16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta

Time-averaged transport in oscillatory squeeze flow of a viscoelastic fluid

Periodically-driven flows are known to generate non-zero, time-averaged fluxes of heat or solute species, due to the interactions of out-of-phase velocity and temperature/concentration fields, respectively. Herein, we investigate such transport (a form of the well-known Taylor–Aris dispersion) in the gap between two parallel plates, one of which oscillates vertically, generating a time-periodic squeeze flow of either a newtonian or Maxwellian fluid. Using the method of multiple time-scale homogenization, the mass/heat balance equation describing transport in this flow is reduced to a one-dimensional advection–diffusion–reaction equation. This result indicates three effective mechanisms in the mass/heat transfer in the system: an effective diffusion that spreads mass/heat along the concentration/temperature gradient, an effective advective flux, and an effective reaction that releases or absorbs mass/heat - in the time-averaged frame. Our results demonstrate that there exist resonant modes under which the velocity peaks when the dimensionless plate oscillation frequency (embodied by the Womersley number, the ratio of the transient inertia to viscous forces) approaches specific values. As a result, transport in this flow is significantly influenced by the dimensionless frequency. On the one hand, the effective, time-averaged dispersion coefficient is always larger than the molecular diffusivity, and is sharply enhanced near resonance. The interaction between fluid elasticity and the oscillatory forcing enhances the efficiency of transport in the system. On the other hand, the identified effective advection and reaction mechanisms may transport mass/heat from regions of high concentration/temperature to those of low concentration/temperature, or vice versa, depending on the value of dimensionless frequency. Ultimately, it is shown that the oscillatory squeeze flow can either enhance or diminish transport, depending on the interplay of these three effective (homogenized) mechanisms

A Parametric Study of Mixing in a Granular Flow a Bi-Axial Spherical Tumbler

We report on a computational parameter space study of mixing protocols for a half-full biaxial spherical granular tumbler. The quality of mixing is quantified via the intensity of segregation (concentration variance) and computed as a function of three system parameters: angles of rotation about each tumbler axis and the flowing layer depth. Only the symmetric case is considered in which the flowing layer depth is the same for each rotation. We also consider the dependence on R, which parametrizes the concentric spheroids (“shells”) that comprise the volume of the tumbler. The intensity of segregation is computed over 100 periods of the mixing protocol for each choice of parameters. Each curve is classified via a time constant, τ, and an asymptotic mixing value, bias. We find that most choices of angles and most shells throughout the tumbler volume mix well, with mixing near the center of the tumbler being consistently faster (small τ) and more complete (small bias). We conclude with examples and discussion of the pathological mixing behaviors of the outliers in the so-called τ-bias scatterplots